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O Abstract 
O 

We study the Laplacian on Stenzel spaces (generalized deformed conifolds), which are 
tangent bundles of spheres endowed with Ricci flat metrics. The (2d — 2)-dimensional Stenzel 
space has SO(d) symmetry and can be embedded in C d through the equation J2i=i z 1 = e " ''■ 
We discuss the Green's function with a source at a point on the S^" 1 zero section of TS d ~ l . 
Its calculation is complicated by mixing between different harmonics with the same SO(rf) 
quantum numbers due to the explicit breaking by the e-deformation of the U(l) symmetry 
that rotates Zi by a phase. A similar mixing affects the spectrum of normal modes of warped 
deformed conifolds that appear in gauge/gravity duality. We solve the mixing problem 
numerically to determine certain bound state spectra in various representations of SO(<i) for 
the d = 4 and d = 5 examples. 



September 2010 



Contents 



1 Introduction [2] 

2 Generalized deformed conifolds [5] 

2.1 Coordinates and metrics [5] 

2.2 Topology and group actions [7J 

2.3 Normalizable functions on Stiefel manifolds M 

2.4 Decoupling of differential equations [10] 

3 The Laplacian on the deformed conifold 1X2] 

3.1 Coordinates and Killing vectors H2] 

3.2 The Laplacian [T4] 

4 Green's Function ITBl 

4.1 Normalization of the Green's function [16] 

4.2 A coupled system of equations [T7] 

4.3 Expectation values of mesonic operators [22] 

5 Basis functions on the deformed conifold |2"B1 

5.1 Group actions on traceless polynomials [2H 

5.2 Homogeneity properties of traceless polynomials [26] 

5.3 Explicit formulae for the Laplacian [23 

5.4 Case studies [28] 

6 Glueball masses for the type IIB background [30] 

6.1 General procedure and case studies [30] 

6.2 Numerical results [32] 

6.3 WKB Analysis [33] 

7 Generalization to Stenzel spaces [351 

7.1 Ricci flat metric and laplacian [361 

7.2 Functions on Stenzel spaces [3H] 

7.3 Green's function computation HQ] 

8 Bound state masses for the CGLP background [43] 

9 Discussion [46] 



1 



A Metric on Stenzel spaces 



1 Introduction 

An interesting class of generalizations of the AdS/CFT duality [ffj3] involves theories which 
exhibit color confinement and discrete spectra of bound states |4j[5] . In such a theory, the 
bound state masses m 8 can be read off from the poles of two-point functions of operators in 
Fourier space, 

(0(k)0(— k)) ~ — r— — — 5- + less singular terms . (1 

^— ' k 2 + mf v / 

In gauge/gravity duality, such poles correspond to normalizable solutions of the linearized 
supergravity equations of motion for the bulk field $ dual to the operator O. 

The first studies of bound state (glueball) spectra via gauge/gravity duality [6]48] were 
carried out for non-supersymmetric backgrounds proposed in R] that describe confinement 
in (3 + 1) and in (2 + 1) dimensions. These backgrounds have SO(5) (SO(6)) symmetry 
due to the presence of an S A (S 5 ) in the geometry. The glueballs that are non-singlet under 
these global symmetries were investigated early on M, with the result that their masses are 
comparable to those of the singlet glueballs. 

The methods of gauge/gravity duality were subsequently applied to calculation of bound 
state spectra in the cascading SU((/c + l)M) x SU(kM) supersymmetric gauge theory coupled 
to bifundamental chiral superfields At, A 2 , Bi, B 2 (5j[To] . This theory has a global SO(4) = 
SU(2)i x SU(2)r flavor symmetry that acts on the chiral superfields. Therefore, it is expected 
that the bound state spectra form Kaluza-Klein towers with increasing SO (4) quantum 
numbers. The gravity dual of the cascading gauge theory is the warped deformed conifold 
(KS) background (5], which is a warped product of 1R 3 ' 1 and the deformed conifold 
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described by the locus J2t=i z i = ^ i n ^ ■ This space has an obvious SO(4) symmetry 
that acts on the ^-coordinates. Therefore, the normal modes of this space can be classified 
according to their SO (4) quantum numbers. The study of glueball masses in the SO(4)- 



invariant sector was initiated in 12 and continued in 13 19]. One of the goals of this paper 
is to initiate a systematic study of glueballs that are not SO (4) singlets (see also |20|). 

We will focus on the case of a minimally coupled scalar field in the bulk, namely a field 
$ that satisfies the 10-d Laplace equation 

□ $ = 0. (1.2) 
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Such an equation of motion arises, for example, from transverse metric fluctuations in the 
IR 3,1 part of the KS geometry. The ten-dimensional type IIB metric of the KS background 



is 



ds 2 = H~^dx ll dx^ + mds%, (1.3) 

where dx^dx 11 is the Minkowski metric on M. 3,1 , ds^ is the Calabi-Yau metric on the deformed 
conifold, and the warp factor H is an SO(4)-invariant function on the deformed conifold. 
Translational invariance in the IR 3 ' 1 directions allows us to look for plane wave solutions to 



eq. (|1.2|) of the form <3> = e lk ' x (f), where <j) depends only on the six internal coordinates. With 



this ansatz, the ten-dimensional Laplace equation takes the form of the eigenvalue problem 

A 6 = - m 2 H(f) , (1.4) 

where m 2 = —k^k^ and A$ is the Laplacian on the deformed conifold. 

When $ is an SO(4) singlet, the dual operator in the boundary theory is the stress-energy 
tensor. The higher SO (4) harmonics of $ correspond to single trace operators which are 
roughly the stress-energy tensor multiplied by polynomials in the bi-fundamental superfields 
of the KS theory and their Hermitian conjugates. The deformation e breaks the U(l)# 
symmetry of the conifold CFT |22| which rotates all the chiral superfields by the same phase 
(in the cascading gauge theory, the U(l)# is broken to Z 2 m by the chiral anomaly |23|, but 
for large M this group may be viewed as a discrete "approximation" to U(l)#). Therefore, 
operators that have the same SO(4) quantum numbers but different U(1)r charges mix in the 
glueball calculation. This infrared mixing of modes in general necessitates solving coupled 
radial equations whose number typically grows with the size of the SO (4) representation. 
This feature adds a new interesting structure to the non-singlet glueball calculations that 
was not encountered in (9l. 

A problem closely related to the Kaluza-Klein spectrum of glueballs is the calculation 
of the Green's function on the deformed conifold. The Green's function is of interest in 
gauge/gravity duality because it determines the extra term in the warp factor produced by 
a stack of p D3-branes on the deformed conifold [24] . The deformed conifold with p added 
mobile D3-branes provides a string dual of the cascading SU((/c + 1)M + p) x SXJ(kM + p) 



gauge theory 25 . As usual, the Green's function may be expanded in an infinite set of 



eigenstates of the Laplacian. We find, though, that the calculation is considerably more 



complicated than in the case of the resolved conifold 26 . For the deformed conifold, finding 
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the eigenstates requires solving the same kind of coupled differential equations that appear 
in the non-singlet glueball calculations, but with different boundary conditions at the tip. 

We will also consider analogous calculations in the weakly curved M-theory background 
found by Cvetic, Gibbons, Lu, and Pope (CGLP) [27]. This background is a warped prod- 
uct of M 2,1 and the SO(5) symmetric Stenzel space 28 , Yl%=i z 1 = f2 ' which is a higher- 
dimensional generalization of the deformed conifold. The CGLP background 27 is similar 
to the KS solution p\, but it is asymptotic to AdS& x without any UV logarithms. In 
the infrared the background contains a blown-up 4-sphere, and the warp factor approaches 



a finite value. Some aspects of the infrared physics were discussed in 29 -31 , but the dual 
infrared gauge theory remains to be elucidated. The CGLP background has a discrete spec- 
trum of normal modes which describe bound states in the dual field theory |32|. They may 
be classified according to their SO (5) quantum numbers, and we will present the bound state 
spectra for a few representations. 

The rest of this paper is organized as follows. In section |2j we explain how one can 



parameterize the Stenzel spaces 28 , Ylt=i z i = ^ \ m terms of a radial coordinate r and 
the coordinates of the undeformed cone with e = 0. The base of this cone is the Stiefel 
manifold V^2 = SO(d)/SO(rf — 2), on which we define a basis of normalizable functions. In 
section |3j we focus on d = 4 corresponding to the 6-dimensional deformed conifold. We 
use the SO(4) symmetry to write the Laplacian Aq explicitly. In section |4j we discuss the 
calculation of the Green's function of this Laplacian with a source placed on the blown-up 
3-sphere. In section [5] we present a complete set of basis functions on the deformed conifold. 
In section |6j we use this basis to derive the eigenvalue problem for the non-singlet glueballs. 
We present our numerical results for the glueball masses belonging to the first few SO (4) 
representations, and compare them with the WKB analysis. In section [7j we generalize our 
results for the Laplacian and Green's function to Stenzel spaces with d > 4. In section [8j 
we compute numerically the glueball masses in the theory dual to the M-theory background 



of 27 . We end the paper with a discussion of our results and of some open problems. 
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2 Generalized deformed conifolds 



2.1 Coordinates and metrics 

The (2d — 2)-dimensional Stenzel space (28) is a Calabi-Yau manifold defined by the quadric 



E 

i=l 



(2.1) 



in C d , where the constant e can be taken to be real and positive. For d = 3 this space is the 



Eguchi-Hanson space 33 , while for d = 4 it is the deformed conifold 11 . We will sometimes 



refer to Stenzel spaces with arbitrary d as a "generalized deformed conifold.' 



When the deformation parameter e vanishes, the space described by eq. (2.1 ) is a cone be- 



cause its definition becomes invariant under rescalings Z{ — > tz-i. The base of this undeformed 
cone is given by 



8=1 



(2.2) 



i=l 



where in order to avoid confusion with the deformed conifold (2.1), we denote the complex 



coordinates in C by yi instead of Zi. The space in (2.2) is also known as the (2d — 3)- 



dimensional Stiefel manifold Vj 2 defined as the set of pairs of orthonormal vectors u and v 



m 



\->\2 \->\ 2 1 

\u\ = \v\ = 1 



u ■ v = 



(2-3) 



These vectors parameterize the coset space SO(d)/SO(d — 2). That (2.2) and (2.3) are the 
same can be seen by the explicit embedding 



Hi = -j= (u>i + ivi) 



(2.4) 



of the Stiefel manifold in C 



The Stenzel space (2.1) can then be parameterized by the yi and a "radial" coordinate 



r: 



(2.5) 
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The constraints (2.2) guarantee that the defining relation (2.1) is obeyed. In order to cover 



the Stenzel space only once, the range of r should be either 0<r<ooor— oo<r<0. 
These two choices are equivalent, so in the rest of this paper we will focus on < r < oo. 
As will be developed more fully in section [2T2 the constant r > surfaces in the Stenzel 
space are much like V^2- they have the same topology as V^2, so on them one can define the 
same functions as on Vd,2- 

As on any Kahler manifold, the metric on this generalized deformed conifold can be 
written in terms of a Kahler potential J 7 as 



3 ' 



9n 



d 2 T 

dzidzj 



(2.6) 



We are only interested in Ricci flat metrics that obey an 0(d) = SO(d) x Z2 symmetry, 
where the SO(d) rotates the Z{ into one anotheiQ and the Z2 acts by flipping the sign of one 
of the Zi. The SO(d) invariance constrains J 7 to be a function of r, so that 



9iS 



e 2 sinh r 



6. 



F"(t)-T'(t) cothr 



e 4 sinh r 



Zi Zj , 



(2.7) 



where the primes denote differentiation with respect to r. The Ricci flatness condition 
reduces to a first order differential equation for J r/ ( r ) given by (27rl 



r 

e 2 sinh r 



d-2 



[21 



The solution can be written as 



d — 2 

with R(t) = — — / (smhv) d ~ 2 dv . 





(2.9) 



Since the Stenzel metric (2.6) has an SO(d) symmetry, it should admit d(d— l)/2 Killing 
vectors. They are 



=T [a) y— +T [a) v- — 



(2.10) 



where the index a runs from 1 to did — l)/2 and the which are real dx d antisymmetric 



1 Because of the SO(d) symmetry (or rather the absence of an SU(c?) symmetry) we make no distinction 
between upper and lower indices on the j/j and Zi coordinates. 

2 For the relation between our conventions and the ones in 27 , see Appendix |XJ See also section [7] below 
for more details on how to obtain (12.8 1. 
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matrices, are the generators of SO(d). One of the reasons that the parameterization (2.5) is 
convenient is that the Killing vectors (2.10) have simple expressions in terms of the y%. In 
these coordinates, the metric (2.7) reads 



ds\ d _ 2 = \T" dr 2 + T' coth r dyidyi 

+ | T' csch r (dyidyi + dyidyt) + ( T" - 7' coth r) yidy^jdyj . (2.11) 

2.2 Topology and group actions 

The meaning of the coordinate r and the relevance of Vd,2 in parameterizing the generalized 
deformed conifolds can be uncovered by examining the topology of these spaces. One can 
start by writing the coordinates z% in terms of the orthonormal vectors U{ and Vi in M. d : 



e (ui cosh — + ivi sinh - j . (2.12) 



If we are interested only in the topology of the generalized deformed conifold parameterized 
by the Zi, we can deform this space smoothly to the set of points in C n given by 

Ui + irvi , (2.13) 

where the range of r is still < r < oo. The real part of this expression parameterizes a unit 
S^ 1 C while the imaginary part parameterizes those vectors in IR d that are tangent to 



it, showing that the Stenzel space (2.1) has the topology of the tangent bundle TS^ -1 ]^ Up 
to reparameterizations, the coordinate r can be identified with the radial coordinate of the 
M. d ~ l fiber of TS . The r = section of the tangent bundle, S' d ~ 1 , is the only non-trivial 
cycle of the generalized deformed conifold. For r > 0, each constant r section is described 
by pairs of orthogonal vectors of fixed length in ¥L d , so it is homeomorphic to V^2- Indeed, 
Vd t 2 is topologically the S d ~ 2 fiber bundle over S d ~ l obtained by restricting the fiber 
in TS^ 1 to the unit S d ~ 2 C R^ 1 . This bundle is trivial if and only if d = 2, 4, or 8 (34||35 
but it has the homology of a product of spheres, S d ~ 2 x S^ 1 , if and only if d is even 36 



We can also understand the topology of a generalized deformed conifold by examining 
the action of SO(d) on this space. The group SO(d) acts on the Stenzel space as isometries 
by rotating the complex coordinates z\ or into each other, while leaving r invariant. This 



3 Thc symplcctic structure induced from C d agrees with the standard symplectic structure on the cotangent 



bundle, T* S d 1 27 . The cotangent bundle T* S d 1 and the tangent bundle TS d 1 are homeomorphic as 



2(d — l)-manifolds, so there is no topological distinction between them. 
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group action is transitive on any constant r surface E T , as any point on that surface can be 
obtained by acting with an SO(cf) matrix O on the point 

z*=e(o ••• isinhf cosh f ) . (2.14) 

When t > 0, the only redundancy in this description is that the matrix O can be multiplied 
on the right by any rotation matrix in the first d — 2 coordinates without changing Oz*. So 
every S T with r > can be identified with the space SO(<f)/SO(<i — 2) = V^2- When r = 0, 
the redundancy is an SO(d — 1) subgroup of SO(d) corresponding to rotations in the first 
d— 1 coordinates, so S has the topology of SO(d)/SO(d — 1) = S^ 1 . 



2.3 Normalizable functions on Stiefel manifolds 

The two problems we are going to address in this paper, namely that of finding the Green's 
function of the Laplacian on the Stenzel space with a source on the S' 1 ' 1 at r = 0, and that 
of finding the eigenfunctions of the Laplacian, involve dealing with functions that are smooth 
at any point on E r with r > 0. The strategy that we will adopt is to use group theory to 
find a basis of normalizable functions on each E r surface, and then expand the functions on 
the Stenzel space in terms of this basis. 

Let us focus on the sections S T with r > 0. While it is certainly true that there is 
some freedom in what the induced metric on S T could beQthe integration measure on this 
space is fixed by the SO(d) symmetry up to an overall r-dependent normalization. In fact, if 
/•tso(<2) is the unit-normalized Haar measure on SO(<i), a canonical way of finding a measure 
/^so(d)/so(d-2) on the coset space SO(d)/SO(d — 2) is through the pushforward operation, 
where each measurable set in SO(d)/SO(d — 2) is assigned a measure equal to the measure 
of its preimage under the projection map SO(<f) — > SO(d)/SO(d — 2). Therefore, up to 
an overall r-dependent normalization, the measure on S T is given by /xso(d)/so(d-2)- As a 
consequence, the space L 2 (E T ) of square integrable functions on S r is the same for all r, 
being equal to the space L 2 {Vd^) with the integration measure Hso(d)/so(d-2)- 



The square integrable functions L (V^) were studied in 37-40 , so most of the results 
presented in the rest of this section are not new. As a vector space acted on by SO(d), L 2 (Vrf i2 ) 
decomposes into irreps of SO(<i). This decomposition can be understood from group theory 
as follows. Any function on V^2 = SO(d)/SO(d — 2) can be thought of as a function on 



4 The most general SO(c?)-invariant metric is a linear combination of Re(dy ■ dy), Im(dy ■ dy), \dy\ 2 , and 
\y ■ dy\ 2 with arbitrary, possibly r-dependent, coefficients. 
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SO(cf) that is constant on each coset, namely f(g) = f(gh), where g is an arbitrary element 
of SO(d) and h is an arbitrary element of SO(d — 2) C SO(d). In order to understand 
the normalizable functions on Vd,2 we should first understand the normalizable functions on 
SO(gT) itself. To construct functions on SO(d), consider an irreducible representation of SO(d) 
on a Hilbert space V with an orthonormal basis e^, 1 < k < dimV. The matrix element 
Mfj(g) = (ei\g\ej) is naturally a function on SO(<f), and by construction it is a polynomial in 
the entries of g whose degree equals the number of boxes in the Young diagram corresponding 
to V. There is a theorem due to Peter and Weyl (see for example [41]) that asserts that 
the set of all such matrix elements between basis states for all irreducible representations 
of a compact Lie group, in particular SO(d), forms a basis for the set of square integrable 
functions on that group, in our case L 2 (SO(d)). Under the left-action of SO(gT) on L 2 (SO(d)) 
(defined as gf(g) = f{g~ l g)), each column in MfAg) transforms like the basis of V, and 
under the right action of SO(d) each row transforms like the dual basis. A consequence of 
the Peter- Weyl theorem is that under the decomposition of L 2 (SO(d)) under the left action of 
SO(d), each irrep V appears with multiplicity dim V, because in the matrix of basis functions 
Mi-{g) there are dim V columns that transform in exactly the same way. 

In constructing a basis of functions for L 2 (Vd.2), we want to retain only those linear 
combinations of MK(g) that are invariant under the right action of SO(d— 2). Since SO(d — 
2) acts only within each row of MX(g), from each row we should keep only those linear 
combinations of basis functions that are invariant under SO(d — 2). It follows that in the 
decomposition of L 2 (Vd j2 ) under the left action of SO(d), each representation V appears a 
number of times equal to the number of singlets in the decomposition of V under SO(d — 2). 

As an example, let us discuss briefly the case d = 4. The five-dimensional Stiefel manifold 



^4,2 = SO(4)/SO(2) is commonly referred to as T 1 ' 1 42 . The group SO(2) representing 
rotations in the first two coordinates of commutes only with rotations in the last two 
coordinates of z*, so it can be chosen to be generated by J^ 3 — Jr 3 , where Jl and Jr are 
the generators of the left and right SU(2) factors in SO(4) = SU(2) L x SU(2) K . For the 
spin jn] irrep, which has a total of (2ji + l)(2j# + 1) basis states, the number of SO(2) 
singlets equals the number of basis states with = mjj, where tul and ttlr are the left and 
right magnetic quantum numbers. Since tul (tur) ranges between dbjx, (±jij), the number 
of such states is 2min(ji, jr) + 1. The discussion in the previous paragraph shows that on 
T 1 ' 1 there must be precisely 2min(ji,j^) + 1 normalizable functions with the same SO(4) 
quantum numbers. Indeed, the fundamental representation of SO (4), [|, |], is represented 
once by yi and once by yf, the [1, 1] irrep appears three times: yiyj, y%\jj-, and y^j — etc. 
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In general, to construct explicitly a basis for L 2 (Vd )2 ), consider polynomials in and y^ 



T\ / rnj2---]n 2 - - 



(2.15) 



7 2 are constant, linearly independent SO(d) tensors that are symmetric in 
the lower and upper indices separately. A choice of basis of square normalizable functions 



where M^../™ 2 



on Vd2 consists of the basis of polynomials in eq. (2.15) modulo the constraints imposed 



by eq. (2.2). In constructing this basis, it is convenient to choose the linearly-independent 



tensors M?^...^ 2 to transform irreducibly under SO(cf). These basis elements can be thought 
of as those functions M^(g) — which, as discussed above, are polynomials in the entries of 
the SO(d) matrix g — that are invariant under the right action of SO(<i — 2). 

As can be seen from the SO (4) example presented above, for a given irrep V of SO(d), 



the linearly-independent polynomials (2.15) with the same SO(d) quantum numbers can be 



distinguished by the total number of yi minus the total number of yi. Later on we will 
call this quantity fh. In terms of the undeformed conifold, fh is identified with half the 
R-charge of the dual operator. Thinking of y\ and yi as independent complex variables, a 



more detailed analysis 38 39 reveals that the polynomials with the same SO(d) quantum 
numbers in V also transform in an irrep of GL(2, C) for which fh is the quantum number of 
a U(l) subgroup of GL(2,C). We will provide more details in section [5j 

To construct wavefunctions 0(r, yi, yA on the deformed conifold, we just take linear com- 



binations of polynomials (2.15) with coefficients dependent on r: 



a 

where F a (y^yi) is a state in a definite SO(d) representation. 



(2-16) 



2.4 Decoupling of differential equations 



The decomposition (2.16) is quite powerful when solving for the Green's function correspond- 



ing to the Laplacian or when solving an eigenvalue problem like in eq. (1.4). The Laplace 



operator on the Stenzel space A2d_2 is SO(d)-invariant, so it can only mix terms in (2.16) 



with the same SO(d) quantum numbers. Therefore, in order to find the coefficient functions 
f a ( T ) f° r a given SO(rf) irrep V, we only need to solve a coupled system of ODEs. The 
number of coupled equations is at most the number of linearly independent functions on 
Va,2 = SO(d)/SO(<i — 2) with the same SO(<f) quantum numbers. As discussed above, the 
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number of such functions equals the number of singlets in the decomposition of V under 
SO(d-2). 

Note that there is no conserved R-charge in the case of the deformed conifold. The would- 
be U(l)i?-symmetry that rotates the z's by a phase is explicitly broken by the non-zero e in 
( 2.1[ ). Nor is there an R-symmetry that rotates the y's by a phase, as can be seen from the 
third term in the metric (2.11 ). Only in the limit of large r, where the e-deformation becomes 
insignificant and the z's and ?/'s are essentially the same, can one define a conserved R-charge. 
Thus, although asymptotically the equations do not couple functions with a different number 
of y's minus y's, there will in general be mixing. Even for finite r however, there is a Z 2 
symmetry that exchanges yi -H- y%. Since it interchanges chiral and anti-chiral operators we 
identify this symmetry with parity. Its existence leads to a decoupling of the functions that 
are even and odd under it. Hence, for a generic SO(d) representation, one has to solve two 
sets of coupled ODEs. 

As a side comment, we note that group theory can be used in similar ways to simplify 
the finding of the Green's function or of the eigenfunctions of the Laplacian on more general 
spaces. The manifold Vd,2 in the above discussion could be replaced by a manifold X on 
which an isometry group G acts with cohomogeneity k (i.e. for which the generic orbits of 
G have codimension k). The space X can be foliated by orbits of points in X under G. If 
K is the stabilizer of a generic point, then the orbit of this point has the topology of the 
space G/K. That G acts on X as an isometry means that each such generic leaf inherits 
from X a measure proportional to the pushforward of the unit-normalized Haar measure on 
G, the coefficient of proportionality being allowed to vary from leaf to leaf as specified by a 
function of the k "transverse" (r-like) coordinates. The space of L 2 -normalizable functions 
on each generic leaf, L 2 (G/K), is an infinite-dimensional vector space acted on by G, so it 
decomposes into irreps of G. As a consequence of the Peter- Weyl theorem, the multiplicity 
of a given irrep V in this decomposition equals the number of singlets in the decomposition 
of V under K. In finding the Green's function or the eigenfunctions of the Laplacian, one 
has to solve in this case a system of PDEs in the k transverse coordinates for each irrep V of 
G. The number of equations in each set equals the number of singlets in the decomposition 
of V under K. Clearly, the number of coupled equations is at most equal the dimension of 
V, and equality holds when each generic leaf in the foliation of X is the group G itself. 
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3 The Laplacian on the deformed conifold 



Before we embark on the discussion of the Laplace operator on the Stenzel space (2.1) with 



arbitrary d, we find it instructive to solve a simpler problem, namely d = 4. There are a few 
advantages to studying this case first. A relatively simple parametrization of this space in 



terms of angles is known 11 , and we will use it to obtain an expression for the Laplacian 



in terms of the yi coordinates; in section [7] we will generalize this expression to arbitrary d. 
In addition, the product structure of the isometry group S0(4) = SU(2)x, x SU(2)# makes 
it easier to understand the group theory that stands behind the decoupling of the various 
harmonics. 



3.1 Coordinates and Killing vectors 

For V42 = T 1,1 , a standard parameterization of the yi coordinates is 
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where 



yi = ^/|( a i & i ~ 

V2 = -j^axbi + a 2 b 2 ) 



&1 i(ix + ± s 

ai = cos — e v 2 4 ' 
bi = cos — e H 2 + 4 ' 



Vz 

2/4 



V2 



(aib 2 + a 2 b\ 



V2 



(oi6 2 - a 2 bi 



CL2 = sin — e v 2 T 4, 

b 2 = sin — e V 2 + 4 ) 
2 



(3.1) 



(3.2) 



The function F'ir) appearing in the metric (2.7) of the deformed conifold is given by 



J-'(r) = e 4/3 (coshrsinhr - r) 1/3 . 



(3.3) 



To describe the six SO(4) = SU(2)^ x SU(2)ij Killing vectors, it is convenient to split the 
generators T^ a ' into two groups corresponding to the left and the right SU(2) factors. The 
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generators of SU(2) L are 



T (i) 



/ 1 o o\ 

-10 

1 
\0 0-10/ 



T (2) 



/ 1 \ 

-1 

-10 

\ 10 / 



T (3) 



and the generators of SU(2) fi are 



T (4) 



/ 10 \ 

-10 

-1 

\ 1 / 



T (5) 



/ 1 o\ 

1 

-10 

\ -10 0/ 



T (6) 



/o -l\ 
0-10 
10 

\1 / 

(3.4) 



/ l\ 
0-10 
10 

y-i ooo/ 

(3.5) 



These generators have been normalized so that they satisfy the algebra [7w, T^] = e l ^ k T^ k \ 
p-(i+3) ;T (j+3)] = e ijfc T (fc+3) ; all Qther commu t a tors being zero. The Killing vectors can be 



computed straightforwardly from (2.10). In terms of the angular coordinates, they are 

9 

I 

9 



£(5) 
^(6) 



902 



9 9 

- COS 2 + COt 6 2 Sill 02 TTT C SC 6» 2 Sin 02 77" , 

O0 2 (702 9?/> 

9 9 9 

-sin 02— - cot6'2cos02 7r— + esc 6 2 cos 2 777 , 

96 , 2 902 9-0 

9 



(3.6) 



9 9 9 

- cos 0i — + cot 9i sin 0i — esc #1 sin X — 

96^1 90i dip 

9 9 9 

- sin 0i — cot 61 cos 0i — — + esc 9i cos 0i — 

9t^i 90i dip 



Quite nicely, the first three Killing vectors act by SU(2) rotations of b\ and b 2 , while the last 
three act by SU(2) rotations of a x and a 2 . We define the quadratic Casimir of SO(4) to be 



(3.7) 



a=l 
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Denoting by 



JLi = i£®, Jm = ^ {i+3) , i = 1,2,3 (3.8) 

the angular momentum operators corresponding to SU(2)l and SU(2)#, respectively, one 
can see that C can be written as 

C = 2(J 2 L + 4) . (3.9) 

From the angular momentum algebra, it follows that functions that transform under the 
b'iji-R] representation of SO(4) have C eigenvalues 

C = 2[j L (j L + l)+j R (j R + l)\ . (3.10) 



We will explain later how to construct such functions. The reason for the factor of 2 in (3.9) 
is that this way the usual SO (4) spherical harmonics, for which ji = j R = £/2, will have 
eigenvalues £(£ + 2). 

3.2 The Laplacian 

The Laplacian on the deformed conifold can be decomposed into a sum of four terms ac- 
cording to their angular dependence: 

A 6 = T + g c (r)C + g n {r)TZ + g c (r)C . (3.11) 

The first term, 

is a second order differential operator in the radial variable r. The others contain derivatives 
only with respect to the angular coordinates. In terms of the complex coordinates jji, they 
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can be written a^3 



d 2 d 2 d 

C = yiVj^~^ \ViVS ~ ^ijVhVk) a a _ + 3y fc - h C.C. 



n 











<9 



9 



1 <9 2 9 

2 (x/»X/j + ViVj ~ tiijUkyk) Qy.Qy. +yk dy^ + C ' C ' 



(3.13) 



The coefficients of these terms are functions of r: 



9c(r) 



2 coth r 



9n{r) 



1 2 coth r 



4 csch t 

r 



(3.14) 



with T 1 given in (3.3). The operator C is the same as in (3.7). It can be checked that C 



1Z, and C each commute with Jl and Jr, so C, TZ, and £ do not mix wavefunctions with 
different SO (4) quantum numbers. 

Our results agree with the formula for the Laplacian on the deformed conifold derived 
up to a typo: we believe that the function A 2 {t) = 2 -10 / 3 coth |(sinh2r — 2-r) 1 / 3 
should be 2~ 7//3 cothr(sinh2r — 2-r) 1 / 3 . The correspondence between our 



m 



24 



defined in 



24 



differential operators C, 71, C, T in equation (3.13) and the operators Or, Di, Dm, and D T 



in 24 is 



C 



-2(Di + n 2 ) -4D 



R 



71 



-4D 



R 



C = D r 



r = a 



(3.15) 



We can also rewrite the coefficients gc, g-R., and gc from (3.14) as 



9c 



1 coth 2 r 
>/ 3 2A{t) 2 ' 



9n 



r4/3 



coth r 



2A(r) 2 AB(t)< 



9c 



coshr 



e 4 / 3 A(r) 2 sinh 2 r 



, (3.16) 



where B(r) = 2" 1 / 3 3" 1 sinh 2 r/(sinh2r - 2r) 2 / 3 , as in 



24 



In writing down eq. (3.131, we do not necessarily require that all coordinates yi and yi should be treated 



as independent — one could solve, for example, for y±, yi, and the real part of j/2 using (2.2 1 before applying 



any of the differential operators in (3.13). But one doesn't have to. Given any two functions fi(yi,y~i) 
and f^itJi^yi) that are equal on the constrained surface (2.2), the formulae in (3.13) are written such that 
T>fi = D/2 on this constrained surface, where T> can be either of the operators C, 1Z, or C. That T>f\ — T>fi 



can be checked by writing T> explicitly in terms of angles using the formulae in section 3.1 and then applying 
the chain rule to convert all derivatives with respect to angles into derivatives with respect to j/j and yi. 
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4 Green's Function 



The results of the previous section enable us to find an expression for the Green's function 
of the Laplacian on the deformed conifold. This question is of interest in gauge/string 
duality because it is relevant to the description of a stack of p D3-branes on the deformed 
conifold, which describes one of the vacua of the SU(p + (k + 1)M) x SU(p + kM) gauge 
theory |5||24)|25). In this case the 10-d metric is again of the form ( |1.3[ ), but the warp factor 
is H = ifxs +SH, where ifxs is the warp factor of the warped deformed conifold solution [5] 
and 5H is proportional to the Green's function with a source at the location (to,?7o) of the 



D3-branes 24 . In particular, if G is the Green's function normalized so that AqG integrates 



to 1 over the deformed conifold, 5H is given by 26 



SH(r,y) = -(2vr) g s p(a') G(r } y;T 0} y ) 



(4.1) 



where g s is the string coupling constant and l/(27ra') is the string tension. For simplicity, 
let us consider the D3-branes placed at a point on the blown up S 3 at r = Op 



4.1 Normalization of the Green's function 



To normalize our Green's function, let us examine the form of the metric close to r = 0. 
With the complex coordinates yi written explicitly in terms of their real and imaginary parts 



as in (2.4), the metric at small r takes the approximate form 



dsi 



2M 
3* 



1 r 2 

-dr 2 + du 2 H dv' 

4 4 



— \- [du 2 + dw 21 \ , 
3s 



(4.2) 



where to obtain the last equality we introduced the coordinate w = rv jl. Geometrically, 
the vector u parameterizes the blown up 5* 3 at r = 0, and for a given u, the vector w 
parameterizes the tangent M 3 plane at u, with \w\ — r/2 being the radial coordinate in this 
tangent plane. 

The equation satisfied by the Green's function of unit-strength source at u = uq and 
r = is 



A 6 G{r, yi ;u ) = — A 5 (3 \u-uo)S (3) (w) . 



(4.3) 



This problem was considered in 24 but our results are different; we find it necessary to include the 



mixing between different wave functions with the same SO(4) quantum numbers. 
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The right-hand side of this equation is normalized so that it integrates to 1 over the whole 
deformed conifold. Note that both u and w were defined as vectors in IR 4 , but each of 
them is forced to live on a three-dimensional surface because of the constraints \u\ 2 = 1 and 



w ■ u — 0. The first delta-function in (4.3) is understood to integrate to 1 over the unit S 3 
parameterized by u with the standard metric on it; the second delta-function is understood 
to integrate to 1 over the IR 3 parameterized by w, again assuming a standard metric. 

4.2 A coupled system of equations 



The solution to the Green's function problem (4.3) can be expanded in terms of polynomials 



in Hi and yi as in (2.15) that have definite transformation properties under SO (4). The 
fact that the delta-function source is at r = as opposed to a generic value of r restricts 
drastically the number of SO(4) irreps that contribute to this expansion. To see which SO(4) 
irreps contribute, one can use the symmetry of the problem as follows. 

The group SO (4) acts transitively on constant r surfaces of the deformed conifold, and 
under the action of SO (4) any point can be brought into the form 

Zi = e (o % sinh \ cosh (4.4) 

for some r. A delta-function source at a generic point (for which r > 0) therefore preserves an 



SO (2) symmetry corresponding to rotations in the first two coordinates in (4.4). This SO (2) 
symmetry is enhanced to SO(3) at r = 0, the SO(3) acting on the first three coordinates in 
(4.4). The Green's function corresponding to a delta-function source at a point on the S 3 



at r = can thus be expanded in harmonics invariant under an SO (3) subgroup of SO (4). 
The embedding of this SO (3) into SO (4) depends on the precise location u of the delta- 



function source, and is generated by the linear combinations of the SO(4) generators (3.8) 
that, acting as Killing vectors, vanish at u = uq. A particularly simple choice corresponds to 
uo = (0, 0, 0, 1), where ( |2.10[ ) implies that the SO(3) subgroup that leaves this point invariant 
is generated by Jl + Jr\\ With this choice, writing Jl and Jr as differential operators as in 
(3.8), one can see that only functions of y 4 and y 4 are invariant under J L + Jr^\ 

For the spin [j'l, representation of SU(2)i x SU(2)^, one can ask how many states 



7 For a source at an arbitrary u — uq on 5* 3 , the unbroken SO(3) is generated by 0( Jl + Jr)O t , where 
O is an SO(4) matrix that rotates (0,0,0, 1) into (uq, Uq, Uq, Uq). The matrix O can be taken to be uil4 



yTV-) — UqT( 2 > — ujT^ 3 \ where I4, is the 4x4 identity matrix and the TW matrices were defined in (3.4). 
^Functions of y± + t/| + Z/3 , ViVi + ViVi + 2/3^3, and y\ + y\ + y| are also invariant, but they can be 



expressed as functions of j/4 and j/4 using the constraints (2.2). 
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are invariant under J L + J R . Clearly, such a state should be invariant under J L3 + J R3 so 
it should be a linear combination of mz,)\jji, — mi). Invariance under Jl± + Jr± fixes 
jh = 3r and all coefficients in the linear combination to be proportional to (— l) mL . So if 
3l 7^ 3r there are no states invariant under J L + J R , and when j L = j R , only the state 
Xlm L =-j L ( — ^) mL \jL, m L)\jL, — m i) is invariant. Quite nicely, this result can be reproduced 
from thinking about polynomials in yi and y~i. As already discussed, invariance under Jl + Jr 
restricts these polynomials to be functions of y& and y& only. It is a straightforward exercise 
to check that Jf = J R when acting on functions Fiy^y^), so all such functions are states in 
some [jl, j R ] irrep with j L = j R . These states are not eigenstates of J^ 3 or J R3 . 



In section 2.3 we explained why for a given spin [jliJr] irrep of SO (4) there are precisely 
2min(jL, + 1 functions on T 1 ' 1 with the same SO(4) quantum numbers. It follows that 
for each SO(4) irrep with = j R = J, there should in fact be (2 J + 1) linearly independent 
polynomials F(y 4 , f/ 4 ) transforming in this irrep that are invariant under + J R . One can 
find them by making the ansatz F{i/ 4) y 4 ) = y 2 4 m or F(y±,y±) = yl m f (y^) , where m 



was introduced at the end of section 2.3 The equation J|F(t/ 4 , |/ 4 ) = 1)^(7/4,^4) then 
reduces to a hypergeometric differential equation for /. Requiring regularity at all points on 
T 1 ' 1 , one finds that the (2J+ 1) linearly independent polynomials are^] 

yf 1 if m > 

F~ 3 ra(y^yi) = Af 3 rh 2-Fi(-J+ \m\, 1 + j+ \m\; 1 + 2\m\;2y 4 y 4 ) x { (4.5) 

y 4 if m < , 

where rh takes values —J, — J + 1, . . . , j, and Aj m is a normalization constant that we take to 
be 

■A/jm = (-l)-J (-2)1*1 V/2J + T f! + |?h • (4.6) 
This normalization is chosen such that 

/ d b y FILFjm = StfSfiri? . (4.7) 



VoKT 1 - 1 ) J T i,i 

This formula is independent of the metric one chooses on T 1,1 as long as the corresponding 



9 



These functions can be expanded in the basis of states described in the following section as 
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volume form is SO(4)-invariant. The hypergeometric function appearing in (4.5) is just a 
polynomial in 2/43/4 of degree j — \rh\; for example: 



Fq,i 



Fi 1 = 2y 4 

2'2 

Fi 1 =2y 4 

2' 2 



F M = 2V3yj 

F lfi = -V3 (1-4^4) 



^V2y\ 

Fa 1 = -4V2y4(l - 3y 4 j7 4 ) 



Fa a 

2 > 2 



Fa_i =-4^^4(1-32/4174) 

Fa _a 

2 ' 2 



(4.f 



4V2yl 



With the normalization in (4.6), we can find how the differential operators (3.13) act on the 
polynomials F^: 



CF- 3 r n = A%j + l)F w , 



jm 

7? F~~ = Arh 2 F~~ 



£Fj„ = (j + m)(j - m + l)F 3 ^i + (j - m)(j + m + l)Fj^+i . 
The Green's function G(r, 2/«; Uq) can be expanded in the basis of F^ functions, 



(4.9) 



G(r,yi;u )= ^ feW^w) , 



j=o ± 1,... m =-J 



(4.10) 



with the coefficients fjm{ T ) being functions of r. Because the source in (4.3) preserves 
the Z 2 symmetry of the deformed conifold, the fjm( T ) are n °t all independent. The Z 2 
symmetry acts on the coordinates by exchanging 2/i 2/i, and therefore maps Fj™(2/4,2/4) — ^ 
Fjm(y4,yi) = Fj_m(2/4, ^4)- Thus, for given J there are rij = l"^^ 1 ] = + 1 independent 
functions fjm( r ) = fj,-m( r ) with rh > 0. 



Using (3.11), (4.7), and (4.9), one can derive a system of coupled ordinary differential 



equations for these functions: 



T fjrh + [4j0 + l)g c + 4m V] fjm 

+0+m + l)(j-m)g c fj,rh+i 

+(j-m + + m)gcfam-\ = V 2 J + 1 



50 



27r 3 e 4 sinh 2 r 



(4.11) 



where the differential operator T is given in (3.12) and gc, gn, an d gc are defined in (3.14). 
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The right hand side of (4.11) follows from 



i d5 y ^ (3) ( M - M o) 5 (3) (^)^(2/4,y 4 ) 



8 S(r) r , 

9e 4 sinh 2 r /V ' 



(4.12) 



and i^(l/V^, 1/V2) = y/2j + 1. In evaluating the integral in the previous equation we 
have assumed that the measure on T 1 ' 1 is normalized so that the volume of T 1,1 has the 
standard value VolfT 1 ' 1 



16^3 

27" 
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. To see the coupling between different components /: 



jin 



explicitly, we write (4.11 ) for j — 0, \, 1: 

Tf = 3 * (T) 
/0 '° 27r3e 4 sinh 2 r ' 

T/i 1 + (3^ c + ^ + 1 



2'2 



3^ 5(r) 
2 2 27r 3 e 4 sinh 2 r 




8fit + 4flfc 25- £ 




3^ 5(r) 
27r 3 e 4 sinh 2 r 



(4.13) 
(4.14) 

(4.15) 



Let us now discuss the boundary conditions one has to impose in solving the equations 
for fjrh- Since for each j we are solving a system of nj = I jj + 1 second order differential 



equations, there are 2nj integration constants. Near r = 0, eq. (4.3) takes the form 



1 d 2 d 

T h 

t 2 dr dr 



-C + TZ + 2C 

2T 2 



G(t, Vi] Mo) 



3 2 / 3 

28/3 e 8/3 



5 (3) (m-%)5 (3) (^) 



(4.16) 



where C, 7?., and >C are n~ 3 x rij matrices. Ignoring the delta-function source, eq. (4.16) has 
power-law solutions that depend on the eigenvalues of the matrix — C + 1Z + 2C. It can be 
checked that these eigenvalues are — A(k — 2)(k — 1) with k = 1, . . . , nj, so there are nj regular 
solutions that behave as r 2 ^ -1 ) and nj solutions that blow up as l/r 2 ^ 1 . Because the first 



term in (4.16) is nothing but a 3-d Laplacian in flat space and the RHS contains a 3-d delta- 



function in w, we infer that G(T,yi;u ) ~ 1/r at small r. The precise normalization will be 
discussed shortly. The requirement that G(T,yi]Ho) ~ 1/r at small r fixes rij—1 integration 
constants to zero, namely the coefficients of l/r 2fc_1 with k > 1. In an asymptotic expansion 



of (4.3) as r — > 00, there are again 2n~ 3 integration constants (of course, related linearly to 



the ones appearing in the r — > expansion), half of them multiplying divergent terms and 
half of them multiplying solutions that vanish asymptotically. Requiring regularity in this 
limit fixes an additional n~ 3 integration constants, namely the coefficients of the divergent 
terms. The last integration constant is fixed by imposing suitable normalization. To sum 
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up, for each j there is a unique solution to the n~ 3 coupled ODEs in (4.3); at small r this 
solution behaves as 1/r with a coefficient that we will now determine. 
The coefficients cjm in the expansion 



fjm{r) 



+ O(t ) 



(4.17) 



must be designed so that the numerical factor and, in particular, the delta-functions on the 



right hand side of (4.16) are obtained when acting with the differential operator of (4.16) 
onto (4.10). The delta-function on IR 3 is readily produced by 



1_^_ 2 d_l 
t 2 dr dr r 



(4.18) 



Since the second term in (4.16) would introduce an unwanted r-dependence, we need to 



choose the coefficients such that the 1/r term of the Green's functions is annihilated by the 
operator — C + 1Z + 2C. It can be checked that the kernel of this operator is spanned by the 
functions 



(4.19) 



m=-j 



Thus by choosing the coefficients Cjm = cj to be independent of m, we eliminate the second 
term and arrive at 



Id „ d -C + K + 2C 



t 2 dr dr 



2r 2 



G{T,yi]u ) 



J=o,|,i,.. 



(4.20) 



In fact, the functions (4.19) depend only on the combination u 4 = -^(3/4 + 2/4), as can be 
seen from the fact that the only normalizable functions annihilated by 

d 2 
-C + K + 2£ = -Av A — + 2(1 -ul- v ^ ' 



dv 4 



Vi, dv 2 



(4.21) 



are independent of W4 = r^(z/4 — 1/4)- If we denote U4 = cos 9, the E~ 3 can be written as 

E 1=V WTT 3[am+ o m . (4.22) 
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For our choice of u = (0, 0, 0, 1), the delta-function on S 3 can be decomposed as 



5 {3) (u-u ) 



7=0,1,1,- 



(2J+1) sin[(2J+l)0] 



2tt 2 



sin 8 



(4.23) 



so the right hand side of (4.20) has just the required structure. Now, we can read off that 



the coefficients Cj„ have to be 



x/2j+T 3 2 / 3 



7T 3 2 8 / 3 e 8 / 3 



(4.24) 



in order to exactly reproduce (4.16). This fully specifies the boundary conditions for the 



Green's function at small r. Note that, upon summing over j, the 1/r term cancels in the 
Green's function except at the location of the source; this is analogous to the cancelation 



of singularities in the Green's function on the resolved conifold 26 . Following the recipe 



described above, it is possible to solve for the functions fjm{r) numerically. Such a numerical 
study is beyond the scope of the present paper. 



4.3 Expectation values of mesonic operators 

Placing p mobile D3-branes on the warped deformed conifold provides the string dual of the 



cascading SU((/c + l)M+p) x SXJ(kM +p) gauge theory 24 25 . The 3p complex coordinates 
of the D3-branes parameterize the mesonic moduli space of the gauge theory; they determine 
the expectation values of the mesonic operators. A special locus on the moduli space, which 
preserves an SO (3) subgroup of the SO (4) global symmetry, is to place all p D3-branes at 
the same point of the S 3 at r = 0. Their effect on the background is encoded in a shift of 
the warp factor proportional to the Green's function on the deformed conifold with a source 
at the location of the stack. The falloffs of various perturbations of the background at large 
r ~ e T l 3 determine the expectation values of various mesonic operators at this point on the 
moduli space. 

From the coupled equations for the radial functions fymif), we find that their leading 
behavior at large r is 

hrn ~ ^r- 2 [ 1+ ^ 1+3J(J+1) - 3 ^ 2/4 J , (4.25) 

where C~ 3 m are constants determined by numerical integration of the coupled ODEs. Since 
the unperturbed warp factor behaves as (In r)r~ 4 at large r, we note that the perturbations 
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of the transverse part of the metric, r 2 H x ' 2 (r)du}^,i t i, fall off as r A J~ A (lnr) 1 I 2 where 



3m 2 

1 + 3J(J+ 1) — -1 



(4.26) 



Up to the power of In r that appears due to the cascading behavior in the UV, this asymptotic 
behavior has the standard form corresponding to a VEV of an operator of dimension Ajm- 



Let us compare eq. (4.26) with the formula for the dimensions of scalar operators dual to 



a mixture of the metric and 4-form perturbations with SU(2)i x SU(2)# x U(l) R quantum 



numbers jl, Jr, R 43 



A 



1 + t^JlUl + 1) + ^JrUr + 1) - 



3R 2 
l6~ 



(4.27) 



The operators described by (4.26) have jz = 3r = J and m = R/2 = —j, — j+ 1, . . . , J. This 



is consistent with the fact that an SO(3) global symmetry is preserved, but U(l)ij is broken 
by the deformation of the conifold. Thus, only the SO(3)-invariant mesonic operators get 



VEVs. The form of these vertex operators in the gauge theory may be obtained from (4.8) 
by replacing y^ with ^(Aii?2 — A2B1) and taking an overall trace. For 



m = j (m = —j) 

these operators are chiral (anti-chiral) with protected dimensions 3|m|, but for other values 
of m the dimensions are not protected by supersymmetry. 



5 Basis functions on the deformed conifold 



In section 2.3, we defined a basis of functions on T ,l that consisted of linearly-independent 
polynomials M^ 2 ..^" 2 y^y^ ■ ■ ■ yi ni y^y^ ■ • ■ yj„ 2 defined up to the equivalence relations im- 



posed by the constraints (2.2). In this section and the next we will find it convenient to 



represent each equivalence class by a manifestly traceless polynomials constructed so that 
the tensor M^ 2 ..^™ 2 is traceless in any pair of the n = rii + n 2 indices. Tracelessness in any 
two lower indices arises because the trace part is proportional to Ylt=i Vi = 0; an< ^ likewise 
for the upper indices. The reason why we impose tracelessness in an upper and a lower index 
as well is that the trace part is proportional to X^=il2/«| 2 = 1> giving a polynomial of lower 
degree. As we will explain briefly, if we restrict ourselves to manifestly traceless polynomials 
we will be able to write down the operator C from section |3.2| in terms of the generators of 
an auxiliary SU(2) group. 



23 



5.1 Group actions on traceless polynomials 



As we have shown in section 2.3, for the \jL,jii) irrep of S0(4) = SU(2)l x SU(2)# there 
are precisely 2min(j£, j#) + 1 linearly independent functions on T 1 ' 1 representing the same 
SO (4) state. Let us try to understand this last statement by examining in more detail the 
traceless polynomials in yi and yi. First, one can check that the operator 



N 



d d 
oyi oyi 



(5.1; 



which counts the degree of a monomial in yi and yi commutes with any of the Ju and Jm 
in the subspace of traceless polynomials. So it is consistent to restrict our attention to a 
subspace of homogeneous traceless polynomials of degree n. 

A traceless tensor M^ 2 ..^™ 2 transforms in general in a reducible representation of SO (4). 
Its irreducible components differ by how many lower indices are symmetrized or anti-symmetrized 
with how many upper indices. We determine the occurring irreducible representations by 
building up M from the product of n\ fundamental representations 4 y , representing the 
lower indices, and n 2 fundamental representations 4y for the upper indices. Writing SO (4) 



as SU(2)l x SU(2)a, one can write 4 % 



1.2 > %\y> 



so the traceless symmetric product of n\ 



4^'s is [^, ^\ y - Similarly, we keep only the traceless symmetric product of n 2 4 s 's, which 
gives the irreducible representation [^, ^h- Lastly, we should only keep the traceless part 
of the product of these representations, which transforms as 



"71 
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n rii — n 2 \ 
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• (5-2) 



A decomposition like (5.2) appears for any partition of the degree n of the polynomials 
(2.15) into two integers, n — ni + n 2 . We therefore conclude that such traceless polynomials 
of degree n are organized into SO (4) irreps [| , | — k] and [| — k, |], where k is an integer 
satisfying < k < l/^l and each representation appears precisely n + 1 — 2k times. In 
particular, C, 71, and £ can mix together up to n + 1 — 2k functions at a time. 

The n + 1 — 2k distinct sets of polynomials differ only in the relative number of y^s and 
f/i's, and we can rotate between them with an SU(2) having the generators 



d d 
oyi oyi 



J-2 



d d 
oyi oy { 



■h 



d 







oyi oyi 



(5.3) 



Here J3 counts (up to normalization) the number of yi minus the number of y^. Its eigenvalue 
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rh running from — | + A; to | — A; is a magnetic quantum number that labels the distinct 
copies of each SO (4) representation. From J\ and J2 we can construct raising and lowering 
operators 

■h = y~, J- = S,^- (5.4) 

Roughly speaking, J + turns one jji into while J_ turns one i/i into Of course, the 
SU(2) action does not correspond to an isometry of the deformed conifold. It is nevertheless 



useful because, given a wavefunction of a particular rh, one can use the operators (5.4) to 
construct those wavefunctions it potentially mixes with. 

One can check that N, j£, J^ 3 , J^, J R3 , and J3 form a maximally commuting set of 
operators in the subspace of traceless polynomials, so they can be simultaneously diagonal- 
ized. The operator J 2 also commutes with all of these operators, but it is not independent 
because 

J 2 = JI + J 2 r-\n{N + 2). (5.5) 

Homogeneous traceless polynomials F a can therefore be described by six quantum numbers 

a = \n; j L , m L ; j R , m R ; rh) , (5.6) 

where n represents the degree of the polynomial, ji and are SU(2)^ quantum numbers, 
j R and rrifi are SU(2)^ quantum numbers, and m is the eigenvalue of J3, or the magnetic 
quantum number of the SU(2). Given such a state, one can construct states with other 
values of and ttir by acting with the raising and lowering operators Jl± = Jli ± iJi2 
and Jr± = Jm ± % Jj®- One can also construct states with other values of rh by acting with 



the J + and J_ operators defined in (5.4). Standard conventions read 



1 

\n;jL,m L ;j R ,m R ;Tfi± 1) = J±\n;j L ,m L ;j R ,m R ;m) , (5.7) 

VJU + 1) -m(m± 1) 



where as can be seen from (5.5) 



71/71 \ 

J(J+ 1) = JlUl + 1)+JrUr + 1) - 2 (2 + V • (5 - 8) 
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For us, the discussion around (5.2) implies 



Tl 

- = max {j L ,j R ) j = vam(j L , j R ) . (5.9) 

5.2 Homogeneity properties of traceless polynomials 

To understand the significance of SU(2), it is useful to describe the homogeneity properties 



of the traceless polynomials on T 1 ' 1 = as discussed in refs. 38 39 . Since we restricted 
ourselves to explicitly traceless polynomials, we can forget about the constraints (2.2), as any 
equivalence class of polynomials with respect to the relations imposed by these constraints 
has an explicitly traceless representative. So we can think of the yi and their complex 
conjugates y~j as unrestricted complex variables. Actually, we can replace y~j by a new 
set of complex variables y"j that are unrelated to yi, and consider all manifestly traceless 
polynomials Pj L ,j R (yi,yj) in the 4 + 4 = 8 independent complex variables yi and yj. When 



both yj = y~j and the constraints (2.2) are satisfied, these polynomials would reduce to the 
functions on the Stiefel manifold V^. 

The representation theory discussed above shows that Pj L ,j R {y%, yj) transforms in an irrep 
of SO (4) whose Young diagram has two rows with ji+jn boxes in the first row and \jl — 3r\ 
boxes in the second row. The reason why only Young diagrams with two rows are possible 
is that any polynomial in or y~j is automatically symmetric, so the only possible anti- 
symmetrizations are between yi and yj, for example yi§2 — y\yi\ one simply cannot write 
down an anti-symmetric product of three or more yi and/or yj. Given such a Young diagram, 
it must be that each of the first \ji — j R \ columns corresponds to either y^ in the first row and 
yj in the second row, or the other way around. Each of the remaining 2min(j^,j^) boxes 
in the first row corresponds to either a y» or a yj, each choice giving a linearly-independent 
polynomial. For this Young diagram there are therefore 2 min(jL, + l linearly-independent 
polynomials with the same SO (4) quantum numbers, in agreement with the discussion from 
section I 



We can now characterize the homogeneity properties of Pj L} j R (y i} yj), following refs. 38 



39 . Thinking of yi and yj as the columns of a 4 x 2 matrix, we consider the GL(2, C) action 

(yi y)^(y'i y') = {y* v)m, (5.10) 

where M is a GL(2,C) matrix. For a given SO(4) irrep, the Pj L ,j R {yi,ijj) transform under 



the action (5.10) as a representation Pj L j R of GL(2, C) because the various Pj L ,j R (yi, yj) have 
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the same total number of y^ and yj. Therefore 

{PslM %)) = (*Wl*, Vj)) P^jniM) , (5.11) 

{Pjhjiiiy'v y'j)j being a column vector representing the components of an arbitrary poly- 
nomial in the [jl, Jr] irrep of SO(4) with respect to a basis of 2min(ji,j/j) + 1 linearly- 
independent polynomials. The group GL(2, C) has a maximally compact U(2) subgroup, 
and Pj Lt j R can be restricted to a representation of U(2) that we also denote by Pj Lt j R . One of 



the main results of refs. 38 39] is that pj L j R is an irreducible U(2) representation with the 
same Young diagram as the [jljOr] irrep of SO(4), namely two rows, one with jh + 3r boxes 
and one with \ ji — Jr\ boxes. As a consistency check, it is easy to verify that the dimension 
of this irrep is 2min(j£, jn) + 1, as expected. 

From the action ( |5. 10 ) with M restricted to be a unitary matrix, it follows that the 
generators of U(2) SU(2) x XJ(1) N are J t for the SU(2) factor and N for the XJ(1) N 



factor defined as in eqs. (5.1) and (5.3), but with the replacement y~i — > yi. The U(2) irrep 
described above restricts to a spin- J irrep of SU(2), where J = min(ji, jn), and to the U(1)jv 
irrep e ld — > e 2temax ^ L '^ R ' , the latter restriction implying that all polynomials in the [jliJr] 
irrep have degree n = 2 max(j£, jn). The values of j and n that we just found from arguments 



about homogeneity properties agree with those found in eq. (5.9) from a more brute-force 
method. 

We should stress that GL(2, C) and its U(2) subgroup do not correspond to symmetries 
of the Stenzel space — one can use the action of these groups on traceless polynomials just as 
a means of keeping track of the homogeneity properties of these polynomials. Indeed, except 
for the U(l) 3 subgroup in SU(2) generated by J 3 and the Z 2 subgroup that exchanges yi 
with fji, the action of U(2) does not even preserve the constraints (2.2). The action of U(l) 3 



is an R-symmetry when e = 0, but this symmetry is broken to Z 2 for the deformed conifold. 



5.3 Explicit formulae for the Laplacian 



With the exception of £ which fails to commute with J 3 , the operators C, 1Z, and £ of (3.13) 
all commute with N, J|, Jl3, Jr, Jr3, and J3. We will therefore fix n, ji, mi, j^, and ttir 
and look at the (2J+ l)-dimensional subspace indexed by m. Restricted to this subspace, 
C will have a fixed eigenvalue, TZ will be diagonal, and £ will generically have off-diagonal 
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entries. In fact the following identities hold when acting on traceless tensors: 



C = 2(J 2 L + J 2 R ) } 11 = 4J| , £ = -2iJ 2 J 3 + JiN . 



(5.12) 



The only non- vanishing matrix elements of these operators are 

(n; j L , m L ;j R , m R ; m\C\n; j L , m L ; j R , m R ; fh) = 2 [j L (j L + 1) + j R (j R + 1)] , 
(n; j L , m L ;j R , m R ] m\K\n; j L , m L ; j R , m R ; fh) = Am 2 , 
{n;j L ,m L ;j R ,m R ;m±l \ £\n; j L ,m L ; j R ,m R ;rh) = Q =F mj a/(JT m)(j± fh + 1) 



(5.13) 



where j = min^,^) as in (5.9). 



5.4 Case studies 

The highest weight state in the j R ] representation of SO(4), which is also a highest weight 
state of the spin J = min(jx,, j R ) representation of SU(2), is 



h 3L,m = (V2 + m) 23 [(V2 + m)(vi + isy A ) - (y 2 + iy 3 )(yi + isy^ 31, 1r] 



(5.14) 



where s = sgn(j L -j R ). It corresponds to F a with a = |2max(j L , j R ); j L , j L ; j R , j R ; j). It is 
easy to check that hj L j R is annihilated by Jl + and J R+ and that it is an eigenfunction of 



C with eigenvalue (3.10). Acting with J_ several times on (5.14) yields the highest-weight 



states in all the copies of the SO(4) irrep [jL,j R ] indexed by fh: 



JL,jR,m 



2J 
j — fh 



x [{V2 + m)(yi + isyt) - (yi + isy A ){y 2 + iy 3 )] 



\0L~3r\ 



(5.15) 



Let us now give a few explicit examples. 



Degree n = 1. In this case, there is only one possible SO(4) representation [|, |]. The 
SU(2) spin is J = 1/2, so there are two sets of states that mix with each other. For example, 
one can consider the highest weight SO(4) states Fi x ^: 



Fx x x =y 2 + iy 3 , Fx x _x = y 2 + iy 3 



(5.16) 
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In this two-dimensional subspace we have 



c = 3, n = i , c = 



1 

1 o. 



(5.17) 



Degree n = 2. In this case, there are three possible SO(4) representations: [1, 1], [1,0], 
and [0, 1]. The representation [1, 1] has SU(2) spin J = 1, so there are three function that 
mix. The SO (4) highest weight states F ljl> ^ fl are 

Fi,i A = (y 2 + iy 3 f , ^1,1,0 = V2(y 2 + iy 3 )(y 2 + iy 3 ) , ^1,1 -i = % + m) 2 • (5-18) 



In this three-dimensional subspace, 



C = 8 , U = diag {4, 0,4}, C 



( V2 ^ 

2V2 2V2 

y y 



(5.19) 



The representations [1,0] and [0, 1] are related by a Z2-symmetry, so we will only focus 
on [1, 0]. Here j = 0, so there is no mixing. The highest weight state is 



^1,0,0 = (2/2 + iy 3 )(yi + iy*) - (2/2 + y 3 )(yi + iy±) 



(5.20) 



Then 



C = 4 , 7£ = , £ = . 



(5.21) 



Degree n = 3. The possible representations are [§,§], [§, |], and [|, |]. In the case [§,§], 
the SO (4) highest weight states F3 3 ^ are 



^§,^1 = (1/2 + «2/3) 3 , ^f,|,i = ^{V2 + iy 3 ) 2 (y2 + iy 3 ) , 
Fs 3 _i = >/3(j/ 2 + iy 3 )(jfc + %) 2 , Fs 3 _s = (y 2 + iy 3 f 



(5.22) 
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The differential operators C, 1Z, and C evaluate to 



C = 15 



K = diag {9, 1, 1, 9} 



C 



( \ 

3V3 4 

4 3^ 

\ ) 



For r3 11 



L2' 2- 



the functions Fz i A corresponding to the highest weight states are 



(5.23) 



Then 



= (V2 + iys) \{V2 + iyz){Vi + zj/4) - (2/2 + il/aXl/i + Z2/4)] , 



(5.24) 



C = 9 



K = 1 



C 



2 
2 



(5.25) 



6 Glueball masses for the type IIB background 



6.1 General procedure and case studies 

As was already discussed, the spectrum of spin-2 glueballs is given by those values of m 2 for 
which there exist normalizable solutions to the ten- dimensional Laplace equation □ $ = 0. 
Under the plane wave ansatz in the non-compact directions x M , this equation reduces to the 



eigenvalue problem (1.4) 



-m 2 H(r)(p 



(6.1) 



where 



H( 



2(g s Ma') 2 f°° , zcothar-l,. , , x1/ „ 
ax „ (smh x cosh x — x) ' 



c8/3 



sinh x 



(6.2) 



and is a function of the internal coordinates r and (yi,yj). The component of </> along a 
state in the spin representation of SO(4) with magnetic quantum numbers and 
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m R has the form 



fm(r)Frn(yi,yi) 



(6.3) 



The equations satisfied by the coefficient functions fmij) are independent of the values of 
m L and m R , but of course depend on j L and j R . In the previous section we derived a matrix 
representation of A 6 in the basis of the highest-weight polynomials Fj L j Rtrh - So the problem 
is reduced to solving a system of coupled ODEs for the coefficient functions /m(T). 

These equations decouple in the UV, where the eigenmodes are approximately those of 
AdSz, x T 1,1 , up to the effects of the cascade. The latter spectrum was identified in earlier 
work 44 45 as follows: a wavefunction i5 carrying quantum numbers |n; jh\ m L\3R\ wr] m), 



behaves asymptotically as 



ci exp 



(A -4)- +...) +^exp 



-A 



T 



(6.4) 



up to powers of r. The normalizable states have q = 0. Then the asymptotics in eq. (6.4) 



correspond to VEV's of operators of dimension 45 



A = 2 + a/4 + Q[j l (j l + 1) + j R (j R + 1)] - 3m 2 . 



(6.5) 



Note that there is no dependence on n,mL,m R . 

The eigenmodes within each SO (4) representation can be decomposed into an even and 



odd sector under the Z 2 parity symmetry that interchanges yi with y^. From (2.5), it follows 



that interchanging yi and yi is equivalent to sending r — > —r. Consequently, the parity of 
the modes is manifested in the parity of the radial wavefunctions fm{T~), where now r is 
allowed to run from — oo to oo. For the [ji, j R ] representation, the number of even modes is 
n e = |~^±I] and the number of odd modes is n Q = I ^p-l , where j = mm(j L , j R ). 

Let us discuss some examples explicitly. For representations of the form [j'l,0], we have 
J = and therefore no mixing. The wave function has only one component, 0(r, yi, yi) = 



fo(T)F jfiiQ (yi,yi), and (6.1) leads to 



Tf + [2j L (j L + l)gc + m 2 H] f = . 



(6.6) 



For the representations [|, „], the wave function has two components, /±i/2, which are cou- 
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pled by (6.1). In terms of ff, 2 = fi/ 2 ± /-1/2 the equations decouple and we have 



Tft 2 + (3^c +9k±9c + m 2 H) /± = . 



(6.7) 



The wave function for representation [1,1] has the three components fx, f$, /_i according 
to j — 1. The first truly coupled equations appear in the even sector between the functions 
ft = fx + f-i and /+ = f : 



T 



ft 



+ 



8g c + Ag n + m 2 H 2y/2g c 



8g c + m 2 H 



ft 



0, 



(6.8) 



while in the odd sector there is only one equation for f\ — fx — f-x' 



Tfc + {8g c + 4g n + m 2 H) /f = . 



(6.9) 



Because of the mixing of modes, in a typical glueball state several operators of different 
dimensions acquire expectation values. For example, for the parity even = jn = I glueballs 
the large r asymptotics, f± ~ e _7T ^ 3 and /g" ~ e - 2 ( 1 +v / 7)/3^ tW^t both operators of 

dimension 7 and dimension 2(1 + \/7) have expectation values. 



6.2 Numerical results 

To solve the coupled equations, it is convenient to use the "determinant method" of |15|. 
Let us focus on the even modes where the corresponding number of equations is n e , the 
discussion for odd modes being entirely analogous. Regularity of the wavefunctions at r = 
fixes half of the 2n e integration constants to zero. We denote the remaining integration 
constants by Sj, 1 < i < n e . For arbitrary Sj, evolving the differential equations to large r 
generically yields non-normalizable wavefunctions. But because these equations are linear, 



the coefficients q of the non-normalizable parts from (6.4) are linearly related to the s,, 
and there exists some n x n matrix 7(m 2 ) for which c% = 'jijSj. This matrix is degenerate 
precisely at those values of m? that solve the eigenvalue problem (6.1 ). 



To normalize the glueball masses, we define a dimensionless quantity m, which agrees 



with the normalization of glueball masses in 15 , by 



95/331/3 
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where H(0) rs 1.13983(g s Ma') 2 /e 8 / 3 . We have computed the first fourteen values of rh of 



each tower in all representations up to [§,§]• The data is available online at 46 . See 



Figure [T] for a summary of our results. As an example of how one reads off the masses of 
various modes, we plotted the determinant of the matrix 7 described above for the even 
modes of the [|, |] representation in Figure |2j 

The masses within each tower exhibit a linear dependence on the excitation number n. 



Furthermore, as was the case in 15 , the leading behavior appears to coincide for all modes. 



This universal behavior at large n has been derived using the WKB approximation 16 . In 
the following section we review their analysis. 



6.3 WKB Analysis 

The Laplacian acting on any eigenfunction <fi can be written as 



d , , n9/ o d G(t) sinh 2 r m 2 H(r) sinh 2 r" 1 
— sinh2r - 2r 2/3 — + v ; - - — + 



0, 



(6-11) 



where G(r) is the angular contribution coming from [gc{r)C + g-ji(r)7l + gc(r)C](j), whose 
exact form will be irrelevant here. The rescaling (j>(r) = (sinh2r — 2r) _1 / 3 ^(r) then puts 
the equation in the Schrodinger-like form 



[d 2 -Q(r)]V = 0, 



(6.12) 



where the effective potential is 

n( , 4 sinh2r 32 
Q{r) 



sinh 4 r [G(r) + m 2 H (r)] sinh 2 r 

3sinh2r - 2r ~ "9"(sinh2r - 2r) 2 ~ 3 • 2 1 /3( s i n h2r - 2r) 2 / 3 



(6.13) 



The Bohr-Sommerfeld quantization condition states that the glueball masses are given by 
solutions to 



a/— Q(x)dx = nil 



7T 



(6.14) 
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SO (4) Representation [j'l, j'r] 



Figure 1: Non-singlet spin-2 glueball masses for the warped deformed conifold [5]. 
Solid/dashed lines represent states with even/odd parity. The representations are ordered 
according to the mass of the lowest state. The numbers being plotted correspond to the 



values of rh as defined in (6.10). 
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Figure 2: The logarithm of |det7(m 2 )| in the range m 2 £ [0,80] for the even sector of the 
[|, |] representation. The glueball ma 
can be seen to organize into 5 towers. 



[|, |] representation. The glueball masses are located at the dips in the plot. The masses 



where n £ Z + and Tj are the classical turning points at which Q(t) = 0. 10 In the limit of 
very large m and upon reinstating numerical factors to match [TBI, we thus derive 



2 . 2 l/3 37r 2 



m ~ 



.£/Yt) sinhr 



(sinh2r-2r) 1 /3 



n 



0.272 n 2 



(6.15) 



7 Generalization to Stenzel spaces 

In section [4] we explained how to reduce the calculation of the Green's function on the 
deformed conifold to solving coupled systems of differential equations. In this section we 



generalize these results to the d > 4 Stenzel spaces (2.1) (for the d = 3 case, i.e. the Eguchi- 



Hanson space [33] , there are further simplifications due to a symmetry enhancement from 
SO(3) to SO(3) x U(l) and the Green's function is known in closed form [47] ) . In particular, 
the d = 5 Stenzel space is relevant to the CGLP solution of the 11-d supergravity equations 



of motion 27 that is dual to a massive vacuum of a (2 + l)-dimensional field theory. 



10 This equation assumes that the potential has sloping walls at both To and ti, which is correct if jl jn- 
If jh = Jr there are modes for which there is no potential barrier at small r. In this case it is standard to 

4 • 



treat tq = as a hard wall, and the RHS of (6.14) becomes nir 
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7.1 Ricci flat metric and laplacian 



The metric for arbitrary d was already given in (2.7). In order to write down the Laplacian 



we need to compute the determinant and the inverse of this metric. However, we cannot 



compute them directly from the expression (2.7), because not all coordinates and 
(i = l,...,d) are independent due to the complex constraint (2.1). Instead, we first have 



to eliminate two coordinates, say z\ and zi, from the metric using this constraint. Letting 
capital indices A run over (a, a) with a = 2, . . . , d, we can write the independent components 
of the metric as 



ds 2 = G A Bdz A dz B , Gab 



G 



ah 



Gab 



where 



G 



ab 



T' z a z b + zizi5 ab T" - T' coth r (z x z a - z a z x ) [z x z h - z b Zi) 



e 2 sinhr 2z\Z\ 



e 4 sinh r 



2zi Z\ 



and G a b = Gba- Now, it is straightforward to compute the determinant 1 ^ 

r 



G = det Gab 



'l) f 



and the inverse metric 



G 



AB 



2e 2 sinh r 



G al 
G~ ah 



2(d-2) 



r 

2z\Z\ 



(7.1) 



(7.2) 



(7.3) 



(7.4) 



where 



G 



ab 



2e 2 sinhr 2 



T" - J^cothr v-^/- 

Za^b o ; ^ / ^ (^i^a Z a ZjJ(ZjZf, Z^Zj 



T'e 2 sinh r 



and G sb = G b ~ a . 

In terms of these functions, the Laplacian is given by 



A 



d 



2d- 2 



Gdzj 



GG 



AB 







dz 



B 



1 9_ 



GG ab — ) + c.c. 
dz 6 > 



d 



(7.5) 



(7.6) 



11 In passing we note that one can now see that the Ricci flatness condition R a i — d a db In \J det G = 



for (7.3) leads to an equation like (2.8) with some constant right hand side. The purely holomorphic and 



anti-holomorphic factors Z\ and Z\ vanish when taking the logarithm and the mixed derivatives. 
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where in this formula, z\ = y/e 2 — z a z a and Z\ = y/e 2 — z a z a are understood to be functions 
of the other variables. In order to write the Laplacian in a manifestly 0(d) invariant form, 
we replace 



d d z a d 

y °L 

dz a dz a zi dzi 



d d z a d 
dz a dz a zi dzi 



(7.7) 



and treat Z\ and Z\ as independent from the other coordinates. The reason for (7.7) is the 



following. We want to act on functions f(zi, Zj) where so far z\ and z\ were dependent on 
za and za- The derivative d/dz a is thus 



df dz x df 

dz a dz a dz x 



(7.i 



and similarly for d/dz a . With (d/dz a )zi = —z a jz\ we find (7.7). Making these replacements 



(7.6) becomes 



A 



2d-2 



4 











' dzj j dzi 



d_ 

'dzj 







3 



dz, 



+ 



r 1 - r coth t 

T'e 2 sinh r 



(7.9) 



d d 

f 1 smlr' r5ij - e 2 cosh r(ziZj + ZiZj) + e 2 (ziZj + ZiZj) ) ^ ^ 



It is obvious that this operator annihilates the constraint (2.1) as every term contains one 



holomorphic and one anti-holomorphic derivative. The fact that (7.9) annihilates the con- 
straints is a prerequisite for this equation to make sense. 



To bring the Laplacian into a form that corresponds to (3.11 ) in the d = 4 case, we pass 



to a "radial" coordinate r and a set of "angles" (yi,yj) that parametrize the base of the 



undeformed conifold. These coordinates are defined in (2.5). Since the y's are even more 



constrained than the z's, it is easiest to work backwards, i.e. to guess an expression for the 
Laplacian in terms of r and y and then to verify the guess by changing coordinates back to 
z. Inspired by the Laplacian Aq on the conifold (d — 4), we find by trial and error: 



*2d-2 



1 j cosh T 

(d-2) cothr<9 r + ^ - -coth 2 rC + -(1 + 2 csch 2 r)ft + 

2 4 sinh r 



7" 

F'-FcothT / 2 _ 1 ^ 1 2 ^ 

^'e 2 sinhr V 2 2 



(7.10) 
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where 



d 2 d 2 d 

„ { d d\(d d 
n = Via I [Via 

1 d 2 d — 2 d 

£ = o (i/t%- + Mi - bijVkyk) q.. Q .. H — — & — + cc. 



(7.11) 



Here C is the Casimir operator of SO(d). Note also that C, 7?., and £ annihilate the constraints 



(2.2) or result in expressions that vanish when the constraints are used. In order to verify 



that (7.10) is equal to (7.9), we need the inverse of (2.5) which is given by 



t = arccosh^j^/e 2 ) 



Vi 



e e ' z. 



a/2 e 2 sinh r 



(7.12) 
(7.13) 



where the r in the second equation is just a shorthand for the right hand side of the first 
equation. We also need expressions for the y- and r-derivatives in terms of the ^-derivatives. 



These are found from (2.5) by the chain rule: 

d e 



dy t a/2 
1 
2 



oz. 



i dzi 
d d \ I 



d T — - cothr ( Zi— + Zi— ] - - cschr ( z { — + z % 



8_ _d_ 

dzj % dzj 



(7.14) 
(7.15) 



where again r is to be substituted by (7.12). This substitution is important when computing 
the second order derivatives. 



Using the functions gdj), giz( T ), and gc(j) defined in (3.14), the Laplacian can be 
brought into the concise form that we had in the d — 4 case: 



A 2d _ 2 = T + g c {r)C + g n {r)n + g c (r)C 



(7.16) 



with 



T = - d {T d - 2 d ) 



(7.17) 



and the other operators as in (7.11). When plugging the d — 4 expression (3.3) for T' into 



T, we recover (3.12). 
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7.2 Functions on Stenzel spaces 



As in the d = 4 case, functions on the Stenzel spaces may be expanded in homogeneous 
polynomials F(yi,y~i) — see eq. (2.16) — with r-dependent coefficients. The tensor 



Mi 



M2—%n 



(7.18) 



in (2.15) now transforms in a representation of SO(d) but otherwise has the same properties 



as before. We want to find its irreducible components. 

We denote SO(d) (vector) representations by their Young tableaux (pi,P2, ■ ■ ■ ,p r ), where 
Pi is the number of boxes in row i and the total number of rows equals the rank r = \d/2\ of 
SO(d). The fundamental representation d is (1,0, ... ,0). The fc-index symmetric-traceless 
representation is (k, 0, . . . , 0). 



Group theoretically the tensor (7.18) is a direct product of the two irreps (n 1; 0, . . . , 0) and 



(ri2, 0, . . . , 0), which represent the two sets of symmetric-traceless lower and upper indices. 



This product decomposes into the following representations 48 ,49 1 (where ri\ > otherwise 
multiply the representations in the other order): 



no 



(n!,0,...,0) ® (n 2 ,0,...,0) 



^2^2(n 1 + n 2 -i-2j,i,Q,...,0) 

j=0 i 



(7.19) 



where i satisfies n\ + n 2 — 2j > 2i > and ri2 — j > i. The representations on the right 
hand side differ by the number of anti-symmetrizations (counted by the index i) and the 
number of traces (counted by the index j) between lower and upper indices. For instance, 
a representation with j = I will contain tensors that are proportional to one delta function 



Mi 



■3V"3n 2 



11— Jn 



8j n . However, we do not want to include any such representation, because the 
corresponding homogeneous polynomial would reduce to one with fewer powers of y and y 
due to the constraint Ylt=i Villi = 1- Therefore, the only representations that we need to 
consider are those with j = 0: 



(ni,0,...,0) ® (n 2 ,0,...,0) 



traceless 



min(ni,n 2 ) 

^2 {ni + n 2 -i,i,0,...,0) 



(7.20) 



8=0 



This sum of representations corresponds to (5.2) in the d = 4 case. 
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Thus, we need to understand representations of the form 



P 

q 



(7.21) 



alias (p, q, 0, . . . , 0). Their dimensions are 

(P - Q + 1)(P + 9 + d - 3)(2p + d - 2){2q + d - 4) 



X 



(d - 2)(d - 3)(d - 4) 2 
p + d - 4\ /g + d - 5 N 
p + 1 A q 



and the eigenvalue of the quadratic Casimir in (7.11) is given by 



C = p{p + d - 2) + q(q + d - 4) . 



(7.22) 



(7.23) 



Just as in the SO (4) case, the representation (7.20) for some given p and q can be realized 



by different choices for n\ and ri2 as long as p + q = ri\ + n^. There are a total of p — q + 1 
copies of the same (p, q, 0, . . . , 0) irrep of SO(d), corresponding to the spin- 2 ^ 2 representation 
of SU(2). The highest weight state of this SU(2) irrep is the polynomial with the largest 
number of y's obtained for n% = p and ri2 = q. In the case d = 4, this polynomial was given 



in eq. (5.14) with m = p = Jl + ji? and n2 = q = \jl — jn\, but it can be generalized to 
arbitrary d by writing 



hp, q = (y 2 + iyzf q [{y2 + m){yi + iy*) - {yi + wa){v2 + iy?)]' 



(7.24) 



Acting repeatedly with J_ we obtain 



p,q,m 



p-q 



p-g 

2 



m 



(V2 + iya) 



(«2 + m 



p-q 

I 2 



x [(y2 + iya)(vi + m) - (yi + iy^iv* + ^3)]' 



(7.25) 



7.3 Green's function computation 

We can now present a generalization of the Green's function computation from section [4] to 
arbitrary d. Suppose we want to solve the Laplace equation with a delta-function source at 
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= 0,...,0), 



A 2d _ 2 G(r, yi ;u ) = 6^(u - ^^(w) . (7.26) 



Here, we used the decomposition y = u + iv as in (2.4) and denoted w = \v as in section |1J 



The right hand side is normalized such that it integrates to unity over the generalized de- 
formed conifold. The Green's function that solves this equation can be expanded in terms 
of polynomials in t/i and y~i with r-dependent coefficients. Since the source preserves an 
SO(d — 1) symmetry that rotates the last d—lyi into each other, the polynomials that ap- 
pear in the expansion of the Green's function also have to be invariant under this SO(d — 1). 
These polynomials are just functions of yi and y\ and consequently they transform in SO(c?) 
irreps with q = 0. For the representation with highest weight (p, 0, 0, . . .), there are p + 1 
such polynomials indexed by the label m ranging between ±p/2. 

One can find the form of the polynomials that appear in the expansion of the Green's 
function by making the ansatz F(yi,y 1 ) = yf m f(yiyi) or F(y 1 ,y 1 ) = yf m f(yiy~i) and solving 
the equation CF = p(p + d — 2)F. The resulting functions are 

/ p p + d \ It/?* ifm>0, 

Fpm(yi,yi) = Nprh 2F1 \m\ - -, — h \m\ - 1, 2\m\ + 1; 2y 1 y 1 ] x 



2 2 J -2H -r ~ ^ n 

' 1 1 if m < . 

(7.27) 



The hypergeometric function in (7.27) reduces to a polynomials of degree f — \m\ in y\y\. If 



we choose the normalization constant to be 



M pm = (-2)N 1 % ■ (7-28) 

(2|m|)! (s - \m\ )! 



we then have 



CF pm =p(p + d- 2)F ? 



pm j 

7? F - = Arh 2 F - 

/v - 1 pm J- pm ; 

£F pih = y--m + lj ^— hm-2 JF Pi ^_i 

+ g + m + l) (^-m-2)F p ^ +1 . 
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The Green's function solving (7.26) can thus be expanded as 



g(t, yi, u Q ) = ^2 Yl fpm( T ) F i*h(yi>yij 



(7.30) 



p=o Th - 



iol some coefficient functions fpm{T~). Using eqs. (3.14), (7.16), (7.17), and (7.29) it is 



straightforward to find the differential equations for the / P m(r). Since the Laplacian on 
the deformed conifold is an SO(cf) singlet operator, it only mixes together functions with 
the same p. The Z 2 symmetry of the deformed conifold and the source further requires 
(r), so for each p we have to solve for n p = I |J + 1 independent functions. 
As in the d — 4 case we will not find these functions explicitly, but instead will give a 
prescription for finding them. 

For each p, there are n p = [|J +1 second order ODEs which give a total of 2n p integration 
constants. Of these, n p are fixed by requiring f prri to vanish at large r. The remaining n p 
integration constants should be fixed by imposing appropriate boundary conditions at r = 0. 



Close to t = 0, the differential equation (7.26) takes the approximate form 







d -C + 1Z + 2C 



T d-2 Q T 



T 



dr 



2r 2 
d- 1 



G(r,yi;u ) 



1 



d - 2 e 2 ( d - 2 ) 



d-2 
d-1 



u 



Uq 



(7.31) 



By arguments similar to those in section [4j at small r we require f prrL to behave as 

1 



fprh 



-d-3 



o 



T 



d-4 



(7.32) 



with the constants c pr a chosen in such a way as to reproduce the delta-function structure 



in (7.31). For a given p, imposing (7.32) requires the c prri to be such that 2^m c pm-^p"i is 



annihilated by the operator — C + TZ + 2C From (7.29), it can be seen that the kernel of this 
operator is spanned by 



E p = 2^ Fprh > (7.33) 

m=-f 

so Cpm = c p is independent of m. More insight into the E p can be obtained by writing 
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-C + K + 2C differential operator 



-C + K + 2C = -2(d - 2) Vl — + 2(1 - u\ - vf)—, 

OVi ov{ 



(7.34) 



when acting on functions of only u\ and V\. The only normalizable functions of U\ and v\ 
that are annihilated by this operator are those that are independent of t>i, so E p indeed only 
depends on U\. Using the explicit formulae for F pih in (7.27), one can write 



E P = T 



d-2\ „d=? 



C p 2 (cos i 



(7.35) 



where u\ = cos 9 and C%(x) are the ultraspherical (Gegenbauer) polynomials. 12 Choosing 



Cpm Cp 



1 2p + d-2 r(d+l) fd-l 1 



327r d - 1 (d-l)(d-3) T(^) V^-2e 2 ( rf - 2 ) 
one can check that the completeness relation 

5 (n - no) = voiT^^^r^^ (cos ^ 



d-2 
d-1 



(7.36) 



(7.37) 



implies that the LHS of (7.31) reproduces the delta-function structure of the RHS of the 
same equation. 



8 Bound state masses for the CGLP background 



The CGLP background 27 , which is an M-theory generalization of the KS background, has 
the metric 



3 2 = H adx^dx 1 * + Hsdsl , 



12 Since E p is only a function of u\ and transforms in (p, 0, . . .) of SO(d), it must be proportional to one of 
the zonal spherical harmonics Z p {6). The zonal spherical harmonics are 



where dtt S d-i is the standard volume element on S' 



d-l 
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where dx^dx^ here denotes the Minkowski metric on M 2,1 , ds\ is the metric on the d = 5 



Stenzel space, and the warp factor is given by 
H(r) 



3 1 m ^ GLP f°° dx 



V2 



sinh x 



(2 + cosha;) 3 / 4 (l + coshx) 5 /2 ' 



(8.2) 



where mccLP is the quantity denoted by m in 50 . As before, we find the glueball masses 
by solving an eigenvalue problem, namely 



-m 2 H(T)(j). 



(8.3) 



The decomposition of A 8 outlined in eqs. (7. 10)— (7. 11 ) reduces the problem to a finite system 



of coupled ODEs. Furthermore, the asymptotic behavior of the eigenmodes in the UV is 



known 43 . The identification of eigenmodes thus proceeds exactly by analogy to Section 6 



To normalize the glueball masses, we defined the dimensionless mass 



m = m- 



"1CGLP 



c3/2 



U) 



We plot the spectrum of glueball masses for the representations (p,q) = (0,0) to (6,6) in 



Figure |3| The numerical values are available at 46 . 

Here too, m 2 exhibits a quadratic dependence on excitation number. To demonstrate 
this dependence analytically, we again write down the generic Laplacian in this geometry, 



d_ 



(2 + cosh r) 3/4 sinh 3 



d 3 7 / 4 sinh 3 \ cosh 3 \ [G{t) + m 2 H(r)] 



2dr 



+ 



1.5) 



where G{r) represents the angular contribution. Then the field redefinition <f>(r) = (2 + 



coshr) 3 / 8 (sinhr/2) 3 / 2, 0(t) puts it in the Schrodinger-like form (6.12) where 



Q(r) 



27 (-7 + 12 cosh r + 3 cosh 2r) coth 2 \ 3 7 / 4 [G(t) + m 2 H(r)} cosh 3 



128 (2 + coshr) 2 

In the large n limit, it can thus be seen that 

Stt' 2 ( r (cosh|) 3 / 2 v^( 



(2 + cosh r) 3 / 4 



(8.6) 



m 2 ~ 



3 7/4 



(2 + coshr) 3 / 8 



n 



m 2 w 3.433n 2 . 



\.7) 
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Figure 3: Glueball masses for spin-2 excitations of the warped M-theory background [27] 
constructed from the 8-d Stenzel space. Solid/dashed lines represent states with even/odd 
parity. The representations are ordered according to the mass of the lowest state. The 



numbers being plotted correspond to rh as defined in (8.4). 
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9 Discussion 



In this paper we have calculated certain Green's functions and normal mode spectra on 
(warped) Stenzel spaces, and discussed the implications of our results for gauge/gravity 
duality. The Green functions determine the back-reaction of an additional stack of branes 
placed on a warped Stenzel space. Their calculation necessitates solving a mixing problem for 
modes that have the same SO(d) quantum numbers but different asymptotic X1(1)r charges. 
Such a problem can only be solved numerically, but we have given a detailed prescription 
for such a calculation, at least in the SO(rf — l)-symmetric case where the D3- or M2-branes 
are placed at the tip of the deformed cone. 

An interesting problem for further work is to replace the branes by anti-branes. This 
provides a description of metastable states in the gauge theories dual to warped Stenzel 
spaces 50,51 . In this context, it is interesting to study the back-reaction of the stack of 



anti-D3 or anti-M2 branes. For anti-D3 branes on the KS background this was considered 
in 



52 , but only when the branes were smeared. Perhaps our calculation of the Green 



function will eventually help attack the more interesting case of localized anti-D3 branes. 

We have also addressed the eigenvalue problems that arise in the calculations of bound 
state (glueball) spectra transforming in the non-singlet representations of SO(cf). So far, we 
have limited ourselves to the spectrum of the minimal scalar equation, which corresponds to 
longitudinally polarized gravitons. It is of obvious interest to extend the glueball calculations 
to other non-singlet modes. 
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A Metric on Stenzel spaces 



In an earlier paper by Cvetic, Gibbons, Lii, and Pope [27], the metric on the deformed 



conifold, (2.11), was written in terms of SO(d)/SO(d — 2) rotation generators. The SO(d) 
rotation generators Lab, A = 1, . . . ,d were denoted by 



'12 



La 



J 2i 



(A.l) 
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where i = 3, . . . , d, and the metric was written as 



ds 2 = \c 2 dr 2 + a 2 a 2 + b 2 ~o 2 + cV (A.2) 



for some specific functions a(r), b(r), and c(r). In order to compare equations (2.11) and 



(A.2), we need the relations 

dy t d yi = \a 2 + \a 2 + v 2 , (A.3) 

i 

|^^| 2 = z/ 2 , (A.4) 

i 

^2(dyidyi + dyidyi) = a 2 - a 2 . (A.5) 

i 

Then it follows that 

a 2 = -(cothr + cschr)^ 7 ' , (A. 6) 

b 2 = -(cothr - cschr)^' , (A.7) 
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